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Abstract 



This work discusses critically the prospects of measuring CP and T violation effects in high energy 
neutrino factories. For this purpose we develop, in the standard framework with three neutrino 
^ ' flavors, simple expressions for the oscillation probabilities in matter that are valid for high E^. All 

I CP violating effects vanish cx E^'^ and are very difficult to detect with high energy neutrinos. A 

. significantly easier task is the determination of the absolute value \5\ of the phase in the v mixing 

' matrix that controls the CP and T violation effects, performing precision measurements of the CP 

CN ! and T conserving part of the oscillation probabilities. 

O 

^ ■ 1 Introduction 

Measurements of atmospheric [||] and solar 0] neutrinos have recently given evidence or 
strong indications that neutrino oscillations exist. These results, together with important 
r-| ! constraints from reactor experiments , give us precious informations about the neutrino 
masses and mixing that we hope will be of great value to develop an understanding on 
physics beyond the standard model. In this work we will assume that the oscillations are 
only between the three known u flavors, and the surprising and potentially extraordinarily 
important results of LSND will be neglected. 

There is currently a very active interest about the planning of future experimental 
studies on u flavor transitions; the possibility to observe CP and T violation effects in 
z/ oscillations (see |jll|, |l3l) is perhaps the most fascinating perspective. Neutrino 
factories have been proposed as a method to provide intense and well controlled beams 
to perform these studies. Two fundamental properties of a neutrino factory experiment 
are the energy of the muon beam, and the neutrino pathlength L. Many proponents 
of a neutrino factory experimental program are advocating high E^^ 0, in fact as high as 
technically possible [E^ ~ 50 GeV or more) and long pathlength [L ~ 3000-7000 Km). 
Other proponents H are advocating a much lower muon energy (i?^ ~ 1 GeV) and a 
shorter pathlength {L ~ 100 Km) (for the possibility of a low energy u factory see 0). 
A critical discussion of the limits and merits of the two options is necessary. 

The study of direct T-violation effects, comparing for example the probabilities for the 
transitions and z/g — z/^, is in principle very attractive [jlO| , however, until u beams 



of extraordinary purity become technically feasible, this study requires the identification 
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of the flavor and electric charge of e^; this is very difScult to do in a very massive detector 
such as those required for these studies. 

The study of CP violation effects suffers because of a fundamental problem: it is essen- 
tially impossible to construct on Earth two CP antisymmetric long baseline experiments, 
because i/'s or z/'s propagate in a medium of electrons and quarks (and not positrons 
and anti-quarks). The effects of the medium on the u flavors transitions are in general 
large, and even in the presence of a CP symmetric fundamental lagrangian one will find 

The fundamental motivation of the "low energy option" is to perform the measurements 
where the asymmetry induced by the matter effects is neglibly small, since the matter 
effects grow with Ej^^this requires low energy neutrinos; however because of the difficulty in 
focusing low energy neutrinos, the smallness of interaction cross sections and the difficulty 
of flavor identifications the experimental challenges are daunting. 

The key point in favour of the choice of a very high for a neutrino factory is that the 
rate of neutrino events, increases oc E^. This impressively rapid growth of the event rate is 
readily understood, as the consequence of two effects: the average energy of the secondary 
neutrinos grows linearly with E^, and to a good approximation ay oc Ey\ moreover the 
angular opening of the neutrino beam shrinks as 7"^ = {E^/m^)~'^ , correspondingly the 
intensity of the v fluence at a far detector increases as oc E'j^. What is not often sufficiently 
stressed is that increasing E^ one has to pay a very high price: the fluence of lower energy 
neutrinos (for a constant number of muon decays) is suppressed oc E~^. For example 
assuming perfect focusing, non polarized muon beam, and approximating ~ 1 the 
fluence of electron (anti-)neutrinos is: 

where is the number of useful muon decays, is the muon mass and G is the step 
function (the fluence vanishes for E^ > E^). Examples of the fluence are shown in fig. |I| 
and ^. The key point is the fact that for Ey much smaller than E^j, the fluence has the 
simple form oc El/E^. 

What is important in the experimental program of course is not the number of v 
events but the size of the effects of the oscillations on the event rates, and actually still 
more important is the size of the new effects that one want to study. The oscillation 
probabilities are suppressed for high E^, and therefore it is not immediately obvious that 
the rapid growth of the event rate with increasing E^ is sufficient to compensate for the 
suppression of the oscillation probability for higher energy neutrinos. One should also take 
into account the fact that larger Ei, means larger matter effects, and therefore requires a 
larger "subtraction" to extract the fundamental CP violation effects from the data. 

The purpose of this paper is to analyse the size of the CP violation effects for high 
energy neutrinos. In this work "high energy" means the energy range where, for a given 
V pathlength L, the transition probabilities decrease monotonically to zero with growing 
E^. This happens for 

\Aml,\L ^ / L N / lA^^sl \ (2) 
2-K V10%my \^3xlO-3eVV 
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that is for E^, larger than few GeV even for the longest possible distances. This is the 
energy range where the proposed high energy neutrino factories will have most of their 
rate. 

As we will discuss in more detail later (see eq. for large the oscillation prob- 
abilities have the following dominant functional dependences on the u energy and the 
pathlength: 

p.,^.^ ~ e;' 

AP,^^^^{CP)^e;' (3) 

AP^^^^^ (matter) ~ E;^ 

The key point is that the CP violation effects vanish rapidly with increasing energy. 
Also important to note is the fact that the fundamental CP violation effects and matter 
effects have the same asymptotic energy dependence, but different dependences on the 
pathlength L. Integrating these probabilities over the expected energy spectrum for a 
neutrino factory far detector one finds the following scaling laws for different signals: 



Rate ~ Ef, L 



3 r-2 

Rate,,^^^^ ~ P^ L° 

ARate,,^.^(CP)~P2L (4) 
ARate^^^;,^ (matter) ~ P° 

The rate of "oscillated events" is approximately independent from the pathlength L and 
grows linearly with P^. The first effect is the result of cancellation between the decrease 
cx of the neutrino fluence, and the growth oc of the oscillation probabilities with 
increasing distance. The P^ dependence is the result of the combination of the decrease 
of the oscillation probability oc P~^, with the growth of the neutrino fluence and cross 
section. The contribution oi CP (or T) violation effects on the event rate is however 
approximately independent of P^ reflecting a cancellation of of the energy depedence 
APy{CP) oc E~^ with the growth of the neutrino fluence and cross section. The CP and 
T violation effects are therefore more and more difficult to observe with increasing P^, 
because the size of the CP violating effects on the rate is constant, while the "background" 
due to the CP conserving part of the oscillation probability increases linearly with P^. 

The simple argument that we have outlined is apparently in conflict with the results 
of previous works that claim W\ |T3 that the largest the P^ the highest the sensitivity to 



the phase 5. The reason for this apparent discrepancy is simple to understand and quite 
instructive. The leading term of the oscillation probability: 

D ^. _ p{0) _ p(0) _ p{0) _ p(0) _ A ^ ^ (r,\ 
J- leading J- u^^Ue ^ lye^f^ V^^Ve ^ Ve^V,, ^lead ^2 ^ 

is equal for all four channels related by a CP or T operation, however the constant 
^lead depends on the value of cos 5. The rate of oscillated events generate by this term is 
independent from L and grows linearly with P^, therefore in principle the higher the muon 
energy the more precisely the constant Aiead can be measured, and cos 5 determined. A 
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significant part of tlie sensitivity (or "reacli" in parameter space) for the phase 6 claimed 
for high energy neutrino factories actually can be understood as simply as the result of a 
a very high precision measurement of the leading term in the oscillation probability. 

There are two important considerations that should however be made. The first one is 
that to transform a measurement of the constant Aiead into a measurement of cos 6 requires 
independent measurements with sufficient precision also of all other parameters in the 
neutrino mass matrix (two squared mass differences and three angles) and a sufficiently 
precise knowledge of the material along the neutrino path. A second consideration, is that 
this measurement has a fundamental ambiguity. All CP violating effects are proportional 
to sin 5. The leading term in the oscillation probability is a CP and T invariant quantity, 
and in fact depends only on the module | sin5|. A measurement of cos 5 that resulted in 
a value different from or 1, would imply the existence of CP and T violations effects 
in the lepton sector, and allow a prediction of their size but not of sign. Of course such 
a result would be of extraordinary importance, however its limitations should be clearly 
understood. 

The paper is organized as follows: in the next section we discuss our conventions for 
the neutrino mixing matrix, section 3 gives a qualitative discussion of the "geometrical 
meaning" of the phase S, section 4 discusses the u mixing in matter, section 5 and 6 discuss 
the u oscillation probabilities in vacuum and in a homogeneous medium. In these sections 
we develop an expression for the oscillation probability in matter as a power series in 
that can be very useful for an understanding of the potential of high energy machines. 
Sections 7 and 8 contain a discussion and some conclusions. Appendix A contains a 
detailed derivation of the most important result of this paper (eq. |32|); additional material 
is in apppendices B and C. 



2 The neutrino mixing matrix 



We will consider in this work oscillations among three neutrinos. The flavor and mass 
eigenstates are related by a unitary mixing matrix U: 



(6) 



For i/'s the mixing is given by the complex conjugate matrix U*. The mixing matrix U 
can be parametrized in terms of three mixing angles {O12, ^13, ^23) and one CP violating 



phase S. We will use the convention suggested in the particle data book |jT6 
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(7) 



where we have used the notation Sjk = sinOjk and Cjk = cos9jk- We need to specify 
a convention for the labeling of the mass eigenstates. We will deflne the state {u^) as 
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the "most isolated" neutrino and as the hghtest between the remaining two states. 
Calhng mi, 1112 and the three mass eigenvalues and defining: 

= ml- m] (8) 

we therefore have that Am^2 is positive by definition, while Am^s can have both signs, 
moreover |Am23| > Am^2- The three mixing angles are then defined in the entire first 
quadrant: Qj], G [0,71/2], while the phase is defined in the interval b G [— 7r,7r]. All 
points in this parameter space represent physically distinct solutions and parametrize an 
experimentally distinguishable "neutrino world". 

For completeness we note that other conventions for the domain of variability of the 
mixing parameters (for the same parametrization of the mixing matrix we are using) 
are possible. Since S13 and b enter the matrix always in the combination S13 e"*'' and 
— Si3 e**^, it is possible for example to enlarge the domain of definition of ^13 to the interval 
[— reducing the interval of definition of the phase: b G [0,7r]; the point (6*13,-151) 
of the conventions used in this paper is then mapped into the point (—^13, \b\). It is 
also common to consider both signs of Am^2 possible. In this case however the angle 
6*12 varies only in the interval [0, |] with the point (6*12, — |Am^2l) of the new convention 
mapped into the point (1 — 6*12, |Am^2l) of our convention. 

In the following discussion it will be sometimes convenient to consider a single quantity 
with the dimension of a squared mass. In these cases we will use the largest squared mass 
difference Am23, as a dimensional quantity and the adimensional ratio 

... = ^ (9) 

Am^3 

to obtain the other squared mass differences: Am^2 = ^12 Am23, Am^3 = (1 + X12) Am23. 
The sign of X\2 is equal to the sign of Am23. 

The Super-Kamiokande data on atmospheric neutrinos indicate that |Am23| is in 
the range 2-5 x 10^'^ eV^ and the angle ^23 is close to ^, while ^13 cannot be large. 
Reactor experiments like Chooz and Palo Verde have obtained stringent upper limits 
on sin^ 26*13, that together with the result of SK tell us that 6*13 is small (sin^0i3 < 0.05). 
The data on solar neutrinos can be interpreted as evidence for oscillations, and give 
information on the the angles ^12 and Amf2 (with a constraint of 6*13, that has to be 
small in agreement with terrestrial experiments). The allowed region in the parameter 
space is composed of discrete regions, only one of which, the so called large mixing angle 
(LMA) solution with Ami2 ~ 10~^-10~^ eV^ and 6*12 close to (but less than) ^, gives us 
a reasonable chance to observe CP violation effects in v oscillations in a standard three 
flavor picture. 



3 The "geometrical" meaning of |(5| 

The three mixing angles have simple "geometrical" meaning in determining the overlaps 
between flavor and matter eigenstates: 

1. The angle ^13 determines how much electron flavor is in the state lz/3), and how much 
is shared between \v\) and \v2) (fractions sin^ 6*13 and cos^ 6*13 respectively). 
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2. The angle ^23 describes how the non-electron content of the state is shared 
between and (fractions sin^ ^23 and cos^^23)- 

3. The angle 612 describes how the electron flavor not taken by the \ is shared between 
the and the \v2) (fractions cos^ ^12 and sin^6'i2). 

It is possible and instructive to consider also the "geometrical" meaning of the absolute 
value of the phase \5\ that enters in the CP conserving part of the oscillation probabil- 
ities. For this purpose it can be useful to use a graphical representation of the mixing 
matrix with "flavor boxes" , this representation has been used before by several authors in 
particular by A. Smirnov [21|. Some examples of this representation are shown in fig. ^. 
Each panel in the figure shows the flavor content of the three neutrino mass eigenstates. 
In the three panels the values of the three mixing angles is identical, but the value of 
the phase 5 changes. The phase \5\ determines how the muon and tau flavor not taken 
by the is shared between the and the \v2)- For 5 = the has the largest 
(smallest) (z/^) component, while for 5 = n the situation is reversed. It is easy to 
see that the K^^^lz^i)!^ overlap grows monotonically with while the |(i^T|i^i)P overlap 
decreases monotonically, and the opposite happens for | (t';i,r|^2) The range of variation 
is determined by the values of the three angles 6^12, ^23 and 6^13. 

In conclusion: the set of overlaps Kz^olz^j)]^ (that is a complete solution for the flavor 
boxes) is equivalent to a perfect determination of the three mixing angles and of the value 
of the phase 5, but with an ambiguity of sign. The mass-flavor overlaps can be determined 
without ever measuring any CP or T violation effects, and therefore the absolute value \S\ 
can be measured without observing any such effect. Of course mathematical consistency 
imply that, if the determination of \S\ differs from the special values or ±7r, then CP 
and T violation effects must exist, and we can predict their existence and their size but 
not their sign. 



3.1 Quasi— bimaximal mixing 

To illustrate our discussion with a concrete example it can be instructive to consider in 
more detail the cases of "bimaximal" and "quasi-bimaximal" mixing. Bimaximal mixing 



corresponds to the values: 9i2 
matrix takes then the form: 



9. 
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0, for the mixing angles. The mixing 





1 

V2 ) 



(10) 



It is well known that in this case the phase 5 is physically irrelevant. In quasi-bimaximal 
mixing we allow for a small non vanishing value of ^13. In first order, that is neglecting 
9\^ and approximating C13 ~ 1, the mixing matrix becomes: 



U 



1 



1 



-i(l + si3 



|(1 - Si3e 

-Kl + ^ise^') ^ 



Si3e 
1 



(11) 
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Note that the values of |f//xi|, \U^2\, \Uti\ and \Ur2\ are not fully determined by the mixing 
angles but can vary in the interval [(1 — si3)/2, (1 + si3)/2]. The extreme values in the 
interval are reached when 5 = or ±7r and the matrix is real. It is interesting to observe 
that "maximum symmetry" (when the four elements |(t'/^,T|i^i,2)| are all equal) is obtained 
when 5 = ±1 and CP and T violating effects are largest. 

For an understanding of the difference between the cases 6 = and 5 = vr it can be 
instructive to look at fig. ^. The figure describes bimaximal and quasi-bimaximal mixing 
when [/ is a real orthogonal matrix. In this case the flavor and mass eigenstates states 
can be represented as two sets of orthonormal vectors in ordinary 3D space. The three 
panels in fig § show the projections in the (z/^,i/t-) plane of the mass eigenvectors. The 
axis is orthogonal to the plane of the paper coming out toward the reader. The left panel 
represents the case of bimaximal mixing: the lies at 45° in the (i/^, Ur) plane while the 
vectors representing i^i and 1^2 are at 45° with respect to the z/g axis, coming out of the 
plane of the figure. The center and right panels show the projections of the vectors when 
6 13 is different from zero. In the center panel has a small component parallel to z/g, that 
is "out" of the plane of the figure (this corresponds to e**^ = 1 or 5 = 0); in the right panel 
1^3 has as small component opposite to the z/g direction, ("into" the plane, corresponding 
to e*^ = — 1 or (5 = vr). Performing a small rotation of the mass eigenvectors from the 
situation of the left panel, one can easily see the effects on the flavor components of the 
other mass states. For 6 = (middle panel) the z/i has (in absolute value) an overlap 
with (ur) larger (smaller) than |, and viceversa for i'2- When 6 = it, (right panel) 
the reverse happens. Allowing the matrix U to be complex, these two discrete solutions 
become the two extreme cases of a continuum of possible different mixings. 

3.2 Boxes and Triangles 

A graphical description of the available information on the neutrinos mixing has been 
introduced by G. Fogli and collaborators in the form of "triangle plots" A "solar 

triangle" plot describes the information about the mass components of |z/e): {|t/eip, |t/e2p, 
|t^e3p}, while an "atmospheric triangle" plot describes the information about the flavor 
components of |z^3): {|?7e3p, l^/^sP, l^^rsP}- Each plot represents a mapping between 
the values of the three u components and the points inside an equilateral triangle, each 
component being proportional to the distance from the point to a side of the triangle. 
The unitarity constraints: 

J2\U^j\^ = l, and Elf^-jf = l (12) 
" j 

are automatically satisfied since from elementary geometry we know that the sum of the 
three distances is a constant. The triangle plots allow to indicate graphically the allowed 
region for the three components. The element |t/e3p is present in both plots and therefore 
there is a consistency check between the analysis of solar and atmospheric experiments. 
The allowed regions in the solar solar and atmospheric triangles carry no in information 
about 6, in the sense that when the allowed regions in both plots shrink to a single point, 
this is equivalent to an infinitely precise determination of the three mixing angles, with 
no information about 6 . 
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More in general one can define six different "triangle" plots, corresponding to the 
three rows and columns of the mixing matrix, that is to the mass (flavor) components 
of each flavor (mass) eigenstate. The set of any choice of three (or more) triangle plots 
(with one case equivalent to the set of three "flavor boxes" ) is sufficient to describe (with 
redundancy) all four mixing parameters including \6\, leaving however ambiguous the 
sign of the phase. To account for the sign of 6 (that is measurable only with the direct 
observation of CP violation effects), one of the triangles must be "doubled". 



4 Neutrino masses and mixing in matter 

The propagation of neutrinos in a medium, differs from the vacuum case. The effects 
of the medium can be taken into account [jl8l considering an effective potential that is 



independent from the u energy. In the study of flavor transitions only the difference 
between the potentials for different flavors is significant, in ordinary (electrically neutral) 
matter one has: 

V = Viu,) - Viu^) = V{Pe) - V{ur) = ^PlGpn, (13) 

where Gp'vs, the Fermi constant and Uf. is the electron density. The effective potential for 
V is the opposite of the v one: V{V) = — l^(z/). For neutrinos traveling close to the Earth 
surface (p ^ 2.8 g cm~^, Ue — 8.4 x 10^^ cm"^) V ~ 1.06 x 10~^^ eV, that corresponds to 
a length 

V-^ ^ 1850 (MJ^J (^) Km. (14) 

{Ye is the number of electrons per nucleon). The effective Hamiltonian for z/'s or T7's 
propagating in matter can then be written: 

niu) = no + n^, n{v) = H;-n„, (is) 

as the sum of the vacuum hamiltonian: 

no = T^U diagK, ml ml] f/t (16) 

and a matter term that in the flavor basis (neglecting a term proportional to the unit 
matrix) has the form: 

iT^m)afS = V 5ae5f3e (17) 

The effective Hamiltonian in matter can be diagonalized to obtain effective squared masses 
values and an effective mixing matrix in matter that will in general be different for u and u. 
The matrices U^, U^, can be parametrized with the form (|^ obtaining the parameters: 
^u'^y ^13'' 5 ^23*^ cind S"^'", and similarly for z7. The solution of this problem involves a 
cubic equation and can be solved analytically [|I^ to obtain the effective parameters as 



a function of the product V Ey] however the solution is sufficiently complex not to be 
particularly illuminating, and is not repeated here. 

One representative example of the dependence of the effective squared masses and 
mixing parameters in matter on the product VEy is shown in figures ^, |^ and (the 
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density p is proportional to the potential V for a constant value of the electron fraction 
Ye = \). Several important features are clearly visible most notably the "resonance" for 
the angle 6^13 at ~ \/S.m2^\/ {2V). More discussion is contained in appendix C. 



5 Oscillation probabilities in vacuum 



The calculation of the oscillation probabilities in vacuum is a well known problem, that 
is briefly outlined here. The evolution equation for a neutrino state is: 



d 

ax 



Ho Ua 



A1 + 



2E, 



U diag[m^, mg, WI3] 



;i8) 



where the Hamiltonian Tio has been written separating a term proportional to the unit 
matrix that can be dropped because it is irrelevant to the transition probabilities, mi, 
777.2 and ms are the u mass eigenvalues, and U is the mixing matrix. The solution of this 
equation for a; = L is: 



z/^(L) = |f/ exp (^-i ^diag [7771,7772,7773]^ [/"^j //^(O). 
The probability for the z/q — > Uf^ transition [a, (3 = e, fi, r) is then: 



(19) 



Y.Upj{e *2idiagK,m2,m^]} f/*. 



JZ^Pj^k^ajUak exp 



L 

^2E. 



(20) 

(21) 



This expression can be expanded more explicitely (here and in all of the following we will 
always consider only non-diagonal transitions a ^ (3, this is no loss of generality, since 
the survival probabilities can be obtained from unitarity) as: 



Al2 a23 

-f [l-cosAi2] + ^ [l-cosA23] + 
±2 J [sin A12 + sin A23 — sin A13] 



413 



[1 -COSA13] + 



where: 



A 



2E,, 



[ml -m])L 



2E,, 



and J is the Jarlskog [2C] parameter: 



J 



J, 



-im[t/ei u;^ u:, U,2] 

Ci3 Si3 S12 C12 S23 C23 sin 6 
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(22) 

(23) 
(24) 

(25) 
(26) 



where we have also given explicitely the expression in terms of the mixing parameters 
used in our convention. The contribution of the first line in equation (^2|) is symmetric 
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under a time reversal (a replacement a ^ (3) or a CP transformation (a replacement 
U ^ U*), while the contribution of the second line changes sign both in case of a T or 
CP transformation (remaining identical for CPT tranformations). The + (— ) in ( p2D is 
valid for neutrinos (anti-neutrinos) when (a,/5) are in cyclic order, that is: (a, /?) = (e, /z), 
(/i,r), (T,e). For {a, (3) in anti-cyclic order the sign must be reversed. The term [sin A12 + 
sin A23— sin A13] can also be rewritten as the product: 4sin(Ai2/2) sin(A23/2) sin(Ai3/2). 



5.1 High energy limit 

It is interesting to consider the oscillation probability in the limit: 



Am^2-^ 



< 1 (27) 



that is in the approximation when the oscillations associated to the longest frequency 
cannot develop. If Ami2 is in the range suggested by the solar neutrino data, the condition 
( p7[ ) will be satisfied in all proposed terrestrial experiments. Developing equation (|22| ) in 
first order in A12 (and using A13 = A23 + A12) one obtains: 

Pi^ya ^ ^p) = ("^23" + ^13") _ ^^^j ^ A£ ^.^ A23 ± 2 J A12 [1 - cos A23] (28) 

In this expression the first term is the dominant one and oscillates with the frequency 
Am23/(2 £^j,); the other two terms are corrections proportional to A12 = Am'l2 L / {2E,^) . 
The last term in ( P5| ) is of great interest because it describes CP and T violations effects. 
These effects oscillate with the same frequency as the leading term, and therefore for the 
detection it is convenient to choose L and E,^ so that \Am2^\L/{2Ei,) = (2n+l) tt with n an 
integer. When this condition is satisfied the CP and T violation effects have a maximum. 
The amplitude of the oscillations of the CP and T violating effects is proportional to the 
Jarlskog parameter, and to Am^g^ therefore the possibility of the detection of CP and T 
violation effects is possible only if three conditions are satisfied: 

1. Amf2 is sufficiently large. 

2. sin 26*12 = 2ci2 S12 is large. 

3. 6*13 is also large. 

The first two conditions are satisfied only if the explanation to the solar neutrino problem 
is the LMA solution. Note also that the amplitude of the CP and T violation effects also 
grows as oc L/E^. 

Examples of the oscillation probabilities P{ve — > v^) and P{Ve ^^t) in the regime 
discussed here can be seen in fig. ^ and in the top panels of fig. § and |10] that illustrate 
the qualitative features discussed above. 



5.2 Very high energy limit 

For very large E^ (keeping L fixed), also the fast oscillations connected with the larger 
|Am23| cannot fully develop, it is then possible to rewrite expression (22) as a power 
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series: 



/ Amno V 



la/3 



Bap ( . 5 I + 



where the constants A^ib and i?Q,/3 are: 

Aap = All + (1 + ^12)' + 

and 



(29) 



(30) 



5,^ = ±8JX12(1 + Xi2) (31) 

Note that Aaf3 is symmetric for CP and T transformations while 5^/3 is anti-symmetric. 



6 Oscillation probabilities in matter 

In this section we will develop some expressions for the oscillation probabilities in matter 
with constant density. The density along the trajectory of a neutrino traveling inside the 
Earth will change slowly, and for the interpretation of real data it will be necessary to 
integrate numerically the flavor evolution equation taking into account these variations, 
however it is a good approximation, sufficient for the purposes of this discussion, to con- 
sider the density constant for all trajectories that do not cross the mantle-core boundary, 
that is all trajectories that have L < 1.06 x 10^ Km. In the approximation of constant 
density the problem of calculating the oscillation probabilities is elementary, in fact one 
can simply use the expressions developed in the previous section with the replacements 
U — > f/m and m| — > Mj, where and the effective mixing matrix and squared 
mass eigenvalues in matter, that can be easily calculated as a function of the parameter 
2 VEi, = 2 \/2GFneEjj. The limit of this approach is that the expressions for the effective 
mixing parameters in matter are complicated, and the results are not transparent. 

In order to gain understanding, have calculated an expression for the oscillation proba- 
bilities that is valid in the limit of large E^, or more rigorously for y = \Am23\L/{AEi,) < 1. 
In this situation it is interesting to write down the oscillation probability as a power series 
in y, generalising equation (p3). The first three terms of this expansion are: 



P(Ua 



= Aai3 

+ Baf3 

+ Cap 

+ ... 



Am^g 


L 








L 


4^. 






L 



AE,, 




(32) 



LV 



mi{LV) 



The quantities Aap-, Bap and Cap are adimensional constants that depend only on the ratio 

X12 = Am?2/^"^23 

symmetry properties: 



Am?2/Am23 and the four neutrino mixing parameters. They have the important 



A 



aP 



+Apa = +A^0 = +Ai^ 



(33) 
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Bal3 
Caf3 



—B/Sa 
+Cf3a 



(34) 
(35) 



Equation (|3^ ) is the main result of this paper, it is derived in appendix A. Note that the 
oscillation probability (in the approximation of high E,^) has been written as the sum of 
three contributions 

1. A leading contribution oc E^"^ that is invariant for CP or T transformation. 

2. AT and CP violating contribution oc E~^ that changes sign both for a time reversal 
and a CP transformation (and is therefore invariant for a CPT transformation) 

3. A third contribution also oc E~^ that is induced by matter effects. This contribution 
is symmetric for a time reversal but changes sign exchanging u with V. 

The adimensional constants A^js, Bajs, and C^^, can be calculated from the neutrino 
masses and mixing: It is remarkable that the first two constants are identical to the 
coefficients in the vacuum case (equations (|30D and (0)). They can be written as: 



(36) 
(37) 

±^Re{{no)lf, [2 (^o)ae (Ho)e/3 - {HIU {Sae + 5/3e)]} (38) 



Bpa = Im[(7^o)Q/3 (^0Ja/3j 



Cj3a 

here TCq is the free Hamiltonian (equation (|T8|) ), e = Am23/(4ii^,^), and the ± sign refers 
to ly (v). It can be checked that adding to the Hamiltonian a term proportional to the 
unit matrix the coefficients do not change. 

The three contributions to the oscillation probability have different dependences on the 
neutrino pathlength L. These dependencies are also simple power laws when L is shorter 
than ~ 2 (in practice when L is shorter than ~ 1500 Km. Developing equation ( p2|) 
for small VL one obtains: 



P(u^ 



A. 



0/3 



+ B, 



af3 



+ c, 




1 



12 
{VLf 



V 



1 - 
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15 



(39) 



The three contributions to the oscillation probability have dependences oc LP' for the 
leading term, oc LP for the CP and T violation effects, and oc for the matter induced 
effects. 

To clarify the simple meaning of this equation let us consider an experiment with a 
fixed baseline L. The oscillation probabilities for the 4 reactions connected by CP and T 
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transformations can be written as a power series in E^^ 



1. 



, A B C 

V V V 

ABC 

U V V 

ABC 

Piv.^v,) = + (42) 

ABC 

Piv.^V^) = + +... (43) 

U V u 

These probabilities have a leading term that is equal for all four channels, and two next 
order terms one (proportional to B) that is due to the fundamental CP violation effects, 
and one (proportional to C) that is the result of the matter effects on the neutrino 
propagation. If one could excavate a tunnel along the neutrino path (to obtain vacuum 
oscillations) the C term in the probability would vanish, while the A and B would be 
modified. For L < 3000 Km the value of A and B are approximately equal in matter and 
in vacuum. 

In an ideal experimental program one could measure all four transitions, and determine 
separately the coefficient A, B and C, however, if we only two channels are experimentally 
accessible (for example z/g and z/g V^) a single (even ideal) experiment cannot 

disentangle the matter effects from CP violations because they have the same functional 
dependence. There are several possible strategies to solve this ambiguity. One solution 
is to perform two experiments with different baselines. The coefficients B and C have 
different dependences on L (for L smaller than 2V~^ the dependences are approximately oc 
Li^ and oc L'^ respectively) , and a comparison of the oscillation rates of the two experiments 
allow in principle to separate the two effects. 

It can appear surprising that for large E^j the oscillation probabilities in matter are 
so similar to the vacuum case, since the mixing parameters differ dramatically from the 
vacuum case. This is the result of some remarkable cancellations. For example it is 
possible to show that the combination of parameters 

si3 C12 S12 C23 S23 sin 5 Am^2 ^"^23 ^"^13 (44) 
is independent from the matter effects, that is 

mat,;/ mat,!/ — •^Vacuum (45) 

This has the important consequence that the CP violating term of the oscillation prob- 
ability is also independent from the matter effects if the u pathlength L is short with 
respect to the three vacuum oscillation lengths {An/lArrtj^]) and the matter length 27r/V. 
In fact: 



APcp = 2J 



sm (^^] sm (^^] sm (^^^ 
\ ^E, I 4E, i \ , 



1 1 

— — J Amjn AmL Aml^ = — — JT (46) 
32 E?, 32 El ^ ^ 
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Similarly the combinations: 

= All i^^u? + (^O' + (47) 

where A-j^^ = —4 Re[Uaj U^j U*i. Upk] are also independent from the matter effects. 

The existence of these cancellations has been discussed in detail in fl^ in the one mass- 
scale approximation (the limit Amfg 0). The general demonstration is technically more 
demanding, however the results can be readily verified numerically. A simpler and much 
more interesting demonstration of these results can be obtained not by the brute force 
approach of comparing the product of explicit expressions the parameters, but using the 
power series expansion outlined in appendix A. 

It can be however interesting to see the "magic" of the cancellation in (^) in an explicit 
example that is worked out in appendix C. 



7 Discussion 



Plots of the oscillation probabilities for the transitions Ve i^fi and Ve — »• 17^ as a function 
of El, for three different values of the pathlength: L = 730, 3000 and 7000 Km in vacuum 
and in matter are shown in fig. ^ and To illustrate more clearly the behaviour 
of the oscillation probability for high E,^, in fig. |TI], [T^ and jTS], we show plots of the 
oscillation probability multiplied by E^. To compute the matter effects we have used 
the exact expression assuming a homogeneous medium with constant density along the ly 
path. Since trajectories with longer L reach deeper inside the Earth where the density is 
larger, the estimated average density is a function of the pathlength: for L = 730, 3000 
and 7000 Km we have used p = 2.84, 3.31 and 4.12 g cm~^, always assuming an electron 
fraction = 0.5. These examples exhibit a number of striking features that we will 
discuss in the following. 



7.1 The measurement of Oi^ 

The focus of this work is on the measurement of CP violation effects, however here are 
included some comments about the measurement of ^13. This measurement is significantly 
easier, and the "optimization" of an experimental program much less ambiguous than for 
the measurement of CP violation effects. It is important to stress the point that the 
"optimum" choice of L and E^ for this measurement will not in general coincide with the 
optimum choice for CP violation studies. 

The key point for this measurement is the existence of a range of E^, where the prob- 
abilities for the z/g z/^ and z/g ^ Uj. are enhanced because of matter effects. This 
enhancement, clearly visible in fig. |^ and |10|, is present for z/'s if Amga > and for z7's 



if Am23 < 0' therefore evidence of a non vanishing 613 should also determine the sign of 
Am23 • '^^^ enhancement of the transition probability is related to the existence of an 



MSW resonance for the angle 613. The effective angle 9^ becomes in fact ^ at an energy 
E ■ 

lAm^l / lAm^l \ /2.8gcm-3\ ^ 
E,,, ~ cos 2^13 -14.1 ^ '".,2 GeV (48) 



2V ■ I 3 X 10-3 eV^ 
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The position of the enhancement of the oscillation however does not coincide with the 
resonance energy but it is at a lower energy -Epeak < Ej-^s- There are three essential points 
about the matter enhancement that should be stressed: 

• The position of the enhancement (that is the value of -Epeak) is determined to a 
good approximation, for small 6*13, only by |Am23| pathlength L. The value 

of -Epeak can bc easily calculated exactly, an approximate formula that describes 
reasonably well the L dependence is: 



E, 



peak 



|Am23| L 
27T + 2VL 



(49) 



note that for small L this coincides with the highest energy where the vacuum os- 
cillation probability has a maximum {E* = |Am23 1-^/27?), while aymptotically (for 

L — > 00) -Epeak ^ Ej-es. 

The size of the enhancement of the oscillation probability, in good approximation 
(for small 6*13) depends only on the pathlength L. In reasonably good approximation 



peak 



vacuum 



VL 



VL + n 



l + -VL + -{VLf + ... 



(50) 



The width of the enhancement region scales linearly with |Am^ 
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Therefore, knowing Am23, given the pathlength L of an experiment, we know a priori 
where (at what E^,) the enhancement will be present, and also how large it will be. The 
value of the probability in the enhancement region is proportional to sin^ 2^13 that is 
unknown, and therefore is not predictable (at least without a successful theory of the u 
mixing) . 



To illustrate these points in fig. |T4| we plot the value of -Epeak (top panel) and the 
enhancement factor (-Ppeak/-Pvac) for Am23 = 3 x 10"'^ eV^ and several values of 6*13 (the 
values of the other parameters are to a good approximation not important). It can be 
seen that when ^13 is small the curves are independent from its value. 

These results can be easily understood qualitatively, since for this purpose it is sufficient 
to approximate Am^2 — 0. In this approximation the oscillation probabilities involving 
electron neutrinos are proportional to the two flavor formulae (see appendix B). The 
oscillation probabilities for the transitions — > z/g are then given by: 



sin^ 6123 



s2 + c2 (1 T E/E,, 



sm 



Am^L 



s2 + 0^1 T E/E,,,) 



(51) 



(where we have used the shorthand notation s = sin 2^13, c = cos 2^13). It is trivial to 
obtain exactly the energy Ei, of the absolute maximum of this probability and its value. 
The position of the maximum in general does not correspond to the the resonant energy, 
when the first factor is largest, because at the resonance the oscillation length becomes 
very long: 

''"vacuum (E res I 



Xr 



sin 2^ 



13 



(52) 
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(where Avac = ^T^Ey/\Am\^\ is the vacuum oscillation length) and the oscillations do 
not have time to develop. The maximum is found at a lower energy, where the mixing 
parameter is smaller but the phase of the oscillations is close to |. 

A detailed analysis of the position of the "peak" in the oscillation probability and the 
value of the probability at the peak shows the existence of a weak dependence on the other 
parameters of the mixing matrix, Am^g; ^12 and 5. Therefore a careful study of the shape 
of the oscillation probability for the enhanced channel (z/g or Ve depending on 

the sign of ^rn^.^ can in principle give information on 5. This line of research is actively 
pursued p^. A fundamental difficulty is that since matter effects are very large close to 



the resonance, they are not easy to subtract, and the uncertainties of the matter density 
profile along the neutrino path can be refiected into effects on the probability of the same 
size as the CP violation effects. 

7.2 The "vacuum mimicking" region 

Looking at fig. ^ one can note a remarkable feature, namely the fact that the oscillation 
probabilities for L = 730 Km and E^, < 0.5 GeV are approximately independent from 
the presence of matter. This phenomenon has been observed before, in particular by 
Minakata and Nunokawa [§, who refer to it as "vacuum mimicking". At first sight 
the closeness of the oscillations in matter and vacuum for energies as large as 0.5 GeV, 
traveling in ordinary matter may seem surprising, since we can expect, and indeed it is the 
case, the effective squared mass values and mixing parameters are modified by the matter 
potential, however these modifications of the oscillation parameters are not refiected in the 
oscillation proabilities if the neutrino pathlength is sufficiently short, namely ifVL<^l. 
In fact it can be proved (see appendix A) that the difference APmatter = -Pmatter — -Pvacuum 
(for a homogeneous medium) can be expressed as a power series in V L and becomes 
negligible for V L small even no matter how different from the vacuum values are the 
effective squared masses and mixing parameters that correspond to the product V E,^. 
parameters. 



This is illustrated in fig. |T5| that describes different regions in the plane (L, E^) where 
the oscillation probabilities for the ly^^r transitions have different qualitative prop- 

erties. The figure is constructed for Am^s = 3 x 10~^ eV^, Am^2 = 7 x 10~^ eV^, and for 
the potential V = 1.06 x 10"'^ eV (that corresponds to the Earth's crust). 

1. The effects of matter on the oscillation probability can be significant only if L is close 
or larger than 2V~^ that is close or to the right of the line labeled a. For growing L 
the matter effects become more and more important (compare fig. |^, |^ and TU) and 
a more and more serious background for CP violation studies. 

2. Curve b is defined by E^, = |Am23| L/(27r), and indicates the position of the high- 
est energy maximum of the transition probabilities, above this line they decrease 
monotonically without further oscillations. 

3. Curve c is defined by the relation E^, = A'ml2 L / {2tt) and gives the highest energy 
where the slow ( "solar" ) vacuum oscillations have a maximum. 
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4. Matter effects are particularly spectacular for Ei, ~ \ Am2^\/{2V) (curve d) when 6*13 
undergoes an MSW resonance. 

5. The matter effects have negliglible effects on the effective squared masses and mixing 
(and therefore no effect on the transition probabilities) if << Am^2/(2^) (curve 
e). For 612 7^ f this corresponds also (adding a factor cos2^^i2) to the location of the 
MSW resonance for the angle 612- 

6. The "high energy region" that is the main focus of this work corresponds to the 
region above the thick dot-dashed lines, when the probability is not oscillating any 
more. 

7. The "vacuum mimicking region" corresponds to the short pathlength (significantly 
shorter than 2V~^ ~ 3700 Km. In order to access CP (T) violation effects where 
they are large, the energy must also be small enough to see the development of 
several oscillations. This qualitatively corresponds to the region delimited by the 
thick dashed line. 

8. Finally the thick curve labeled A indicate the energy where the MSW enhancement 
(or suppression) of the probability is most important (correspond to -Epeak discussed 
in the previous subsection. The line is plotted only when the enhancement of the 
probability is larger than 20%. The enhancement becomes more and more important 
with growing L (see equation (0). The best place to search for a non vanishing ^13 
is close to this line. 



7.3 High Energy Neutrinos 

We will now discuss in more detail the behaviour of the oscillation probabilities for high 
E,y, when they are monotonically decreasing. Some important features of the oscillation 
probability in vacuum are illustrated in fig. [Tl| thas shows the product Pv^^u^, x El plotted 
as a function of E^^ for a fixed value L = 730 Km: 

1. For large E^, the oscillation probability is well approximated with the form A/E^. 

2. The value of the constant A (keeping all other parameters fixed) depends on the 
value of the cos 5 (this is a crucial remark). 

3. The CP violation effects (present when sin 5 7^ 0) have an energy dependence oc E~^. 

Fig. |T2| illustrates how the oscillation probability is modified by the matter effects (for 
the same L = 730 Km and the same u masses and mixing as the previous figure). The 
probabilities in vacuum and in matter for the longer pathlength L = 3000 Km are shown 
in fig. Some important points are the following: 

1. Matter effects generate a u/v asymmetry. 

2. The matter induced asymmetry has the same energy dependence (oc E^J^) as the one 
generated by the fundamental CP violation effects. 
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3. The effects of the phase S and of matter can both contribute to the observable u/u 
asymmetry, either adding or subtracting from each other. 

4. The effects of matter depend strongly on the pathlength L. 

5. The relative importance of the matter induced asymmetry with respect to the asym- 
metry generated by the fundamental CP violation effects grows approximately lin- 
early with L (see equations (0) and (|39D ), and the "background" of the matter 
induced asymmetry is a much more serious problem at L = 3000 Km. 

6. The effect of the presence of matter on the leading order term {A/E'^) of the transition 
probability also depend on the pathlength L 

— For short L the constant in matter A^^at is equal to the vacuum value A^^c- 

— For longer L [L > 2 V~^) the leading term of the oscillation probability is sup- 
pressed (see again equations (^2])); however one has Amat = F{LV) x A^ac, that 
is the constant A is proportional to the vacuum value with a proportionality 
factor that depends only on the product of V L; therefore a measurement of the 
leading term again carries information about cos 5. 

7.4 The leading term in the oscillation probability and \5\ 

The leading term in the oscillation probability ~ A^^/El is equal for all the four tran- 
sitions related by a CP or a T transformation, however in principle its measurement, 
together with a precise determination of the other oscillation parameters, can give infor- 
mation about the phase 6. In fact most of the sensitivity to 5 claimed by recent analysis 
of the potential of high energy //-factories, is essentially the result of a high precision 
measurement of the constant A. Considering all other parameters as fixed and ^rn^^ pos- 
itive, the constant A^.^ is largest when 5 = 0, and decreases monotonically with growing 
reaching a minimum value for \5\ = tt. Conversely A^t is maximum for |5| = vr and 
minimum for 6 = 0. For Am^s < the dependence of and A^r on \6\ is reversed. 

It is simple to give a qualitative explanation for this behaviour. The constant Aajs can 
be expressed as the sum of three contributions each one associated to a squared mass 
difference: 

A^p = Ail X + X (1 + ^12)' + X 1 (53) 

In this expression each factor A-j^^ = —4:Re[Uaj U*}. U^j U^k] is weighted proportionally to 
to the square of the relative squared mass difference Am^;^. The general expressions for 
the different terms in our parametrization of the mixing matrix are easily calculated . For 
the Ve ^ i^fi transitons one has: 

^11 = 4 S12 C12 (C^3 - 5^3) + (0^2 - 5^2) Si3 S23 C23 COS S (54) 
^11 = 4Si3 C?3 S23 C12 + 4 S12 C12 Sl3 Ci3 S23 C23 COS 6 (55) 
^ej = 4 Si3 Ci3 S23 - 4 Si2 C12 Si3 0^3 S23 C23 COS 5 (56) 

One can see that the term Al^ is non vanishing also for S13 = 0, and in fact is related 
to the "solar oscillations". Note also that the sum Al^ + A^^ is independent from 8 
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(and from 612), however the two individual terms do depend on the phase (and also on 
612). For cos 5 —* 1 the term Al^ becomes largest while A^^ becomes smallest; this is a 
consequence of the fact that the overlap is largest for cos 5 = 1 (see discussion in 

section ^). This has simple but important consequences when considering the combination 
(|5B[). When X12 is positive (that is for Amgs positive) one has jAm^gl > |Am23| and the 
contribution A^^ receives the largest weight, therefore the e-/i oscillations have the highest 
probability when Al^ is largest that is when cos 5 = 1. For X12 negative {Aml^ < 0) one 
has lAmfgl > |Amf2l and therefore it is the term A"^^ that receives the larger weight, 
therefore the probability is largest when ^423 is largest, that is for cos 5 = 0. A similar 
discussion can be performed for e-r transitions. In this case the oscillation probability is 
largest when cos 5 = for Am23 > 0, and when cos 5 = 1 for Am23 < 0. 

For a simple illustration we can consider the case of quasi-bimaximal mixing (^23 = 
^13 = f )• Keeping only terms in first and second order in S13 one obtains: 

Ae^^,er = ^ " ^Is) + 2 S^g + S13 COS 6 [{1 + Xu? " 1] (57) 

One can recognize three main contributions that can be considered as the "solar contri- 
bution" , the "^13 contribution" and the "mass splitting effect" . 

1. The first term due to the effect of oscillations involving only the states ui 
and 1^2 ■ This contribution is non- vanishing also for ^13 = 0, it is "garanteed" to exist, 
since it is responsible for the oscillations of solar neutrinos. 

2. The second contribution oc 5^3, can be understood as oscillations between the state 
Us and the quasi-degenerate pairs z/i-z/2. This contribution is independent from the 
solution of the solar neutrino problem and can in fact exist also for Am^2 0) 
however it can be arbitrarily small since it is proportional to 5^3 for which exists 
only an upper limit. 

3. The third contribution is oc S13 0:12 cos(5, and is non-vanishing only if both Am^2) 
and Si3 are non zero. It arises as the consequence of the different weights for the 
oscillations "between" the pairs of states 1^1-2/3, or 2/2-1^3 due to the associated squared 
mass difference. This is the contribution that carries information about the phase S. 

An important remark is that, for a large interval of values, when S13 decreases it 
becomes easier to measure the contribution of the cos 6 term, even if it becomes smaller. 
This happens because the leading contribution is usually the "6^13 term" that is oc sf^ and 
decreases quadratically with S13, while the cos 5 contribution decreases only linearly with 
S13. In fact Cervera et al in their analysis of the sensitivity of //-factory experiment 



have found (see fig. 22 in their work), that the range of Am^2 where it is possible to 
distinguish 6 = and 5 = 7r/2 is apoproximately constant (with actually a small increase) 
when ^13 becomes smaller, down to the smallest angles they investigated. This observation 
seems a paradox, since the CP violation effects decrease linearly with S13. However the 
observation is correct, and can be easily understood with the argument outlined above. 



If the authors of [0 had studied the possibility to discriminate between S = and 



(5 = +1 the allowed region would be significantly smaller, and shrink linearly with 613. 
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This in fact reveals the conceptual difference between measuring CP violation effects and 
measuring the phase 5. 
In summary: 

1. The measurement of the leading term in the oscillation that is a quantity symmetric 
under CP and T transformations gives information about the value of cos 5. This 
result does not allow to determine the sign of the CP violation effects. 

2. A precise measurement of the leading term in the oscillation probability is significant 
as a measurement oi \5\ only if the other oscillations parameters can be measured 
separately with a sufficient accuracy. The requirement on the precision of the mea- 
surements of the "solar parameters" 9i2 and Am^2 particularly stringent (they 
can be deduced from eq. 53-57). 



8 Conclusions 

In this work we have addressed the question of the optimum strategy, that is the best 
choice of u energy and pathlength, to measure the two remaining completely unkwown 
parameters in the neutrino mixing matrix [^], that is the angle and the phase 5. 

For the measurement of 6'i3 it is possible to make a strong case for a high E^^ and 
long L program. The oscillation probability for the z/g <-> (or z/g depending on 

the sign of Am?) transitions will be enhanced in a well defined and precisely predictable 
energy range. It is in this energy range that the search for the effects of a non vanishing 
^13 has the best possibilities. The maximum of the oscillation probability corresponds to 
an energy E^^ak that is proportional to |Am23|, is only weakly dependent on Am^2 ^^^d 
the mixing parameters, and grows with increasing L in a well defined way (see eq. HUI and 
the following discussion). For lAm^gl = 3 x 10"^ eV^ -Epeak is approximately 1.6, 5.2 and 
7.4 GeV for L = 730, 3000 and 7000 Km. The value of the probability at -Epeak can be 
predicted as P^'^^^ ~ sin^ 2^13 sin^ ^23 x F where F is a matter enhancement effect that 
to a good approximation depend only on the pathlength L (see eq. ^). For L = 730, 
3000 and 7000 Km the enhancement F is ~ 1.25, 2.8 and 10.3. The best strategy for a 
detection of ^13 is therefore to use a very long pathlength (since it improves the signal 
to background) and design a neutrino beam with maximum intensity in the energy range 
where the probability is predicted to have the maximum. Note that for this study, a 
conventional beam ^6, ^ could be competitive with a neutrino factory. 



The study of CP and T violation effects is a fascinating subject, and it is remarkable 
that if two conditions are satisfied: (i) the LMA solution is the explanation of the solar 
neutrino problem, and (ii) ^13 is sufficiently large, these effects are in principle observable 
with accellator u beams, and the phase 6 is experimentally measurable; however these are 
extraordinarily difficult tasks, and the best strategy is not easily determined. Very likely 
in this case the very well controlled and intense beams of a z/ factory are a uniquely well 
suited tool, however the choice of F^ and L is not a simple decision. A large E^ allows 
very high event rates, but results in high energy neutrinos that have small oscillation 
probabilities and for which the CP violation effects are strongly suppressed (oc E~^), 
moreover the matter effects become a more dangerous source of background. 
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To analyse quantitatively if the high rates of a high energy neutrino factory are suffi- 
cient to extract information about the phase 6 we have analysed in detail the oscillation 
probabilities for high energy neutrinos, expressing the probability as a power series expan- 
sion in the adimensional parameter y = |Am23| L/{AEy). This expansion is useful when 
y is less than unity, that is for larger than few GeV, even for the largest possible L, 
and reveals several important features of the probability. Keeping only the lowest order 
terms in the expansion, the oscillation probabilities is the sum of three contributions: 

r2 ^3 ^4 

P = Ao (cos 5)—± BcP sin 5 — ± C^,, V — (58) 

V V V 

One distinguish: 

1. a leading order contribution oc -Ej"^, that is symmetric under CP and T transforma- 
tion 

2. a CP and T antisymmetric contribution of order that depends linearly on sin 5, 

3. a matter induced contribution, also proportional to E~'^ ^ that is invariant for a T 
reversal transformation, but changes sign replacing v with V (or viceversa). This 
term is proportional to the potential V and vanishes in vacuum. 

It is important to note the L dependence of the three terms: longer L enhances the CP 
violation effects, but enhances more dramatically the matter effects. 

The largest effect of the phase 5 is on the leading term coefficient A. This coefficient 
can be measured with great precision, and for this purpose a very high energy is the 
optimum solution, however the value of cos 5 extracted from the measurement only if the 
squared mass differences and mixing angles are determined (from other measurements) 
with sufficient precision. The possibility to obtain the required accuracy for the "solar" 
parameters Q\2 and Am^2 is problematic and should be critically analysed. 

The detection of a genuine CP violation effect, is more difficult. It requires first of 
all to subtract the matter induced asymmetry. This can be done having two experiments 
with different baselines, and using the different L behaviour of the two contributions, or 
studying the energy dependence of the probability down to lower Ey. Note that the effect 
on the event rate induced by the fundamental CP violation is constant with increasing 
since the increase in the rate oc is compensated by the suppression in the probability 
oc E~'^ for higher energy neutrinos. Since the backgrounds increase with energy, the 
optimum solution for the search of the asymmetry is not an arbitrary high energy. 

The authors of some recent works [0, |T3[ on the sensitivity of high energy neutrino 
factories for the determination of 5 have neglected to consider the entire interval of def- 
inition of the phase 6 G [— vr, vr] , studying only the positive semi-interval. It would be 
interesting to see a reanalysis of those works that considers the entire interval of definition 
for b. The outcome should be that at least in a significant part of the parameter space 
an input value ^mput can be reconstructed with fits of approximately the same quality as 
^rec ~ ^input and (5rec ~ " f^input- This would rcvcal how much of the sensitivity is coming 
from the measurement of the CP violating part of the oscillation probability, and how 
much is coming from the high precision measurement of the CP conserving part. 
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An ambiguity of sign in the dermination of 5 leaves ambiguous also the sign of all CP of 
T violation effects in the neutrino sector. This can obviously be a hmitation, for example 

for a discussion of the relation between these effects and the observed CP asymmetry of 
the present universe. More in general, it should be noted that a measurement of cos 6 is 
not rigorosly speaking a measurement of CP violations at all, but it implies the existence 
of CP violations of predicted size but unknwon sign. One point that in my opinion 
would require more attention is the following, the measurement of cos 5 (different from 
the special values and l)is equivalent to the statement that three flavors, two squared 
mass differences and three mixing angles are insufficient to describe all observed results 
about the i/ flavor transitions, and that the inclusion of a new parameter can reconcile all 
results. It is not clear if this interpretation of the data would be unique. It is interesting 
to discuss if other mechanisms, for example the introduction of new neutrino properties 
(such as FCNC interactions), or additional small mixings ( "LSND-like" ) with light sterile 
states, could also be viable descriptions of the data. 

An experimental program with low energy neutrinos and a short pathlength has in 
priciple some very attractive features: (i) the problem of disentangling the matter effects 
is much less severe because these effects are small, (ii) a direct measurement of CP 
violation effects is possible, (iii) the oscillation probability can have more structure, and 
the CP violation effects can be very large (with APcp/P ~ 1). Unfortunately, the 
experimental difficulties are enormous, because of poor focusing, small cross sections, and 
the experimental difficulty of flavor determination. The question of which one of the two 
options (high energy or low energy) is more promising for the measurement of the phase 
5 and of CP violation effects for leptons remains in my view still open, and more detailed 
studies have still to be performed. 
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Appendix A: Oscillation probabilities as power series 



A.l Oscillations in vacuum 

In this appendix we will show how the v oscillation probabilities both in vacuum and in 
homogeneous matter can be expressed as power series. 

We can start with the vacuum case. In this discussed in section ^ the oscillation 

probabilities can be calculated with simple analytical formulae and are a function of the 
ratio LjE^. It is useful to define the adimensional quantity y: 

It is straightforward to expand the oscillation probability as a power series in y: 

oo 

P^^f,{L,E;) = Y.cr^ y- (60) 

n=2 

developing in a Taylor series the trigonometric functions in the analytic expression (^). 
The constants c°^^ are adimensional quantities that can be written as functions of the 
squared mass ratio X12 = Amf2/Am23 and of the 4 mixing parameters. The expansion 
(|60[) is formally always valid, it is of course useful only when y is less than unity, that 
is for short z/-pathlength or high Ey. Since the expansion of the cos (sin) function that 
describe the CP conserving (violating) part of the probability has only even (odd) powers 
of their argument, the coefficients have the following symmetry properties: 

' c^^/^ = +c^^" = +c°^^ = +c^^" for n even, 

(61) 

cT^P = -cf^" = -c"^/^ = +cf^" for n odd, 

^ IL It lb lb I 

Note in particular that the lowest order (leading) term of the expansion is exactly CP 
conserving, while the next term is CP anti-symmetric. It follows that the determination 
in an experiment with a fixed baseline, that the fiavor transition probability in vacuum 
vanishes at high energy with the form Py^^y^ ~ 0'E~'^ + bE~^ with b ^ 0, would be a 
proof of the existence of CP violations in the neutrino sector. 

It is useful to rederive the expansion ( |60D of the oscillation probability with a different 
method, that can be more easily extended to the case of neutrinos propagating in matter. 
The S matrix for fiavor transition can be calculated as: 

S{ua Z//3) = SfBa = exp[-iHo L]po^. (62) 

where Tio is the free Hamiltonian. The S matrix for the transitions of z7's can be obtained 
replacing Hq with the complex conjugate Hq. Expanding the exponential one has: 

Sfia = exp[-iHoL]po^ = 5pa + {-iL)(Ho)i3a + -^{-iLYi'HDpa + ■■■ (63) 
The transition probability can be obtained squaring the corresponding matrix element: 

Po,^f3 = \S(Se,\^ (64) 
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Collecting all terms proportional to for each integer n, one can then obtain the prob- 
ability as a power series in L. This actually corresponds to the powers series in since 
L always enters in the combination Ti^ L or Ti^ L, and the Hamiltonian can be written 
(neglecting a term proportional to the unit matrix) as: 

no = ^U diag[-(l + 2xu). -1, l]U^ = ^h (65) 

and each power of Hq (or 7-^q) contributes a factor Am|3/(4£',^). Writing explicitely the 
lowest order terms one finds: 

Pp^^ = [V Kp] y' + Mkp Ch%^] y' + --- (66) 

from where we can read the expressions for C2^^ and c^~^^ . The symmetry properties of 
the coefficients can be easily checked: 

1. The lowest order term (oc y"^) of the probability is symmetric for time reversal, since 
7^0 is an Hermitean matrix: (7io)a/3 = O^oYpa- 

2. The leading term is also symmetric for a. CP transformation since in this case one 
has to replace the Hamiltonian Tio with its complex conjugate. 

3. The next to leading term oc is antisymmetric under a CP or T transformation as 
can be immediately deduced from the fact that Tip is hermitean. 

The coefficients j/ia/jp and \m\h*^p {h^) ^p] can be easily written in terms of the mixing 
parameters and the squared mass ratio Xi2 verifying that the expansion ( p^ ) is identical 
to the one given in equation 



A. 2 Oscillations in matter 

The effective Hamiltonians describing the flavor evolutions of z/'s and I7's in matter can 
be written as: 

HyL = HoL + peV L = hy + pez (67) 
HuL = n*QL-peVL = h*y-peZ (68) 

where we have introduced the Ue projection operator pe (with the property (pe)" = Pe) 
that in the flavor basis has the components {pe)a/3 = ^ae^ise and a second adimensional 
quantity: 

z = VL (69) 

Writing the S matrix in an expanded form, squaring the element S/sa and collecting all 
terms proportional to (y" z™) one can obtain the transition probability as a power series 
in both y and z: 

oo oo 

P^^piL, Ey) = P^^piy, z) = J2i: 2/" (70) 

n=2 m=0 
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Naively one could expect to see terms of order (y° z"^), [y z) and [y z^) but they are present 
only on the diagonal of the S matrix and are irrelevant for the transition probabilities. 

Note that the set of coefficients c"^^ are of course identical to the vacuum expansion. 
From this we can deduce the very important fact: if z is small, that is when the v 
pathlength is much shorter than the matter length V^^, the oscillation probabilities are 
approximately equal to the vacuum case. This is a sense is not an entirely obvious fact, 
because a condition on the pathlength does not say anything about the importance of the 
matter effects on the neutrino masses and mixing. It is indeed possible that -Ejy and V are 
such that the mixing parameters are entirely different from the vacuum case (for example 
one could sit on a MSW resonance), however, if L <^ the oscillation probabilities 
will coincide with the vacuum case. This result, especially in the 2>i' case, appear as 
the consequence of some remarkable "cancellations" between the effective values of the 
mixing parameters and squared masses, however in this formalism it is entirely natural 
and obvious. 

Some important symmetry of the coefficients are the following: 



n,m 




1 „a^/3 ^ 




for 


n 


even and m even 


„a-^/3 _ 
n,m 




_fM^'P - 
n,m 


n,m 


for 


n 


even and m odd 


„a^/3 _ 
n,m 


_„/3^a _ 
n,m 


_„a^l3 _ 
n,m 


' n,m 


for 


11 


odd and m even 


r<^-*P - 


_^P-^a _ 




n,m 


for 


n 


odd and m odd 



(71) 



Note how the matter effects when they enter with an odd power of the potential have 
opposite signs for u and V. 

An explicit calculation of the coefficients of lower order in y gives for n = 2: 

C2,mfodd) = {^'^) 



while for m even one has: 



CWcvcn) - [^"/3 ^a/j] ^2,m (73) 



"2,m(cvcn) L"'Q!P ''a/JJ "'2,m 

where (i2,m are simple numerical coefficients: 

m (-1)'' 2 

^2,m(even) — ^ 2^ 71 , i M / T^i TTT ~ 7 , n\t ^ \^^) 

^qIk + ly. [m + 1 — k)\ (m + 2)! 
For the elements with n — 3 one has: 

<Xen) = (h'U] d3,n^ (75) 



1 

6 

where the quantities ds^m are numerical coefficients: 



Cs.mfodd) = R [2 hae Kp - {h^)aP (Sae + 5/3e)]} (^S.m (76) 



m (-1)'' 2 

[k + ly. [m + 2 — ky. (m + 2)! 
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da m = 12 1-^-' y TTT = ^'^^'^\^} ^""^ for m odd (78) 

^o(A; + l)!(m + 2-A;)! (m + 3)! ^ ^ 

(The coefficients have been defined so that 1^2,0 = c?3,o = c^3,i = !)• 

It is useful, for example when considering an experiment with a fixed baseline, to 
resum over all z terms, and express the probability again as a power scries in y (that 
corresponds then to a power series in E~^)^ with coefficients that are distance dependent: 
For the lowest terms the sum is easily obtained: 

00 00 /2\2 /' 7\ 

Y.d2,mZ-^ E d^,mZ-'=(-^ (79) 



m=0 m(even)=0 

48 



E ^2,m --^ 



g)-|sin(.) 



(80) 



m(odd)=l 

The reason to keep separate the sums of the m-even and m-odd contributions, is because 
they have different symmetry properties under a CP transformation. For m even (odd) 
the c'^~^ coefficients have the same (opposite) sign for z/'s and z7's. 

It can be interesting to verify the results derived above in the special case of two 
neutrino mixing, when simple exact expressions for oscillations probabilities exist. This 
is done in the next section. 
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Appendix B: Two flavor mixing case 



It can be useful to to discuss v^. or z/g ^ oscillations in the case of two flavor 

oscillations when the explicit calculation of the oscillation probabilities in matter is very 
simple. Moreover this case is an exact solution in two interesting limiting cases: 

1. The limit x\2 — >■ (Am^2 small). In this case one has to replace Q — >■ ^13, Am^ — >■ 
Amig. The oscillation probabilities for the z^e ^ ^y. and v^. ^ transitions are then 
given by the two flavor formula multiplying by sin^ ^23 and cos^ ^23- 



2. The limit ^13 — > 0. In this case one has to make the replacements ^ — > ^12, Am^ — > 
Am^2- The oscillation probabilities for the v^, ^ and v^. <-> v^. transitions are then 
given by the two flavor formula multiplying by cos^ ^23 and sin^ ^23- 



In the two flavor approximation the vacuum oscillation probability is: 

"Am^L" 



-Pvac(i^e i^n) = sin 2^ sin 



To have the oscillations in matter we simply have perform the replacements 6 
Am^ (Am^)in, with: 



(81) 
9m and 



sin^ 29m = sin^ 29 



sin^ 29 + cos 29 ^ 



2E^V' 
Am? 
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and 



(Am^)^ = Am^ \/sin2 2^ + (cos 2^ ^ 



2E^V 



(82) 



(83) 



The minus (plus) sign refers to neutrinos (anti-neutrinos). It can be seen that in general: 
(i) the oscillation probabilities in vacuum and in matter can be very different from each 
other; (ii) the oscillation probabilities in matter for u and I7's are also in general very 
different. 

We are interested in the probability for large E^,. For vacuum oscillations it is straight- 
forward, to expand the probability as a power series in y = Arri^ L / [AE^,): 



P ~ 

^ vac — 



sin^ 29 



' Am^i y sm^29 f Ajr^V 



+ ... 



(84) 



To compute the same expansion for the probabilities in matter it is useful to introduce 
the variable: 

using the shortened notation s = sin 26*, c = cos 26*, the oscillation probability in matter 
can then be rewritten as: 



mat 



^2ce + e' 



sm 



VL I 

V1t2c£ + £2 

2 



(86) 
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Developing in a power series in e and writing for simplicity a 



^ one obtains: 



that can be rewritten as: 

+ 
+ 



s e sin a ± 2 see (sin a — a sin a cos a) + 0{e 



sin^ 26 



(87) 





[1 .21 




— sm a 


[ 4i5. j 


.a 



sin2 2^cos2^ f Am^ V 
3 [ AE, ^ 



sm a — a sm a cos a 



We can note that the first (second) term is symmetric (anti-symmetric) for the replace- 
ment V — i> —V that is replacing v with V. 

Developing the expression (|88|) for small a (that is for L < V~^) and reinserting the 
definition one obtains: 

2 r 



Pn 



sin^ 26 



'Am'V 
AE,, 



1 



+ 
+ 



sin2 2^ cos 2^ ( Ani^ V 



12 



VL 



{VLf 
15 



+ 



(89) 



We have obtained with an explicit calculation a set of interesting results for the oscillation 
probability in the limit of large energy 

(i) The oscillation probabilities for z/'s and z/'s become asymptotically equal, and vanish 
with a leading contribution of order y'^ ~ 

(ii) The matter effects generate a correction of opposite sign for z/'s and F's. This cor- 
rection vanishes more rapidly with increasing energy oc y'^ ~ E'J^ 

(iii) For small pathlength: L < the matter effects grow linearly with the matter 
potential V . 

(iv) Again for small pathlength the leading order term of the oscillation probability grows 
with the pathlength oc L^, while the correction due to the matter effects grows more 
rapidly oc L^. 

Comparing equations (|88D and ( pOD with equations (|32D and ( p9D one can see that we 
have reproduced with an explicit calculation the results for the leading order term of the 
oscillation probability and for the matter effects. It is also easy to see that in the two 
flavor case considered here the free Hamiltonian, and its square can be written in the 
flavor basis as: 

Aw? 



Tit] 



AE, 



-cos2e sin 2^ 
sin 20 cos 20 





" 1 


' 







1 



(90) 



Calculating the coefficient Ae^^, -Be^t and CepL according to the general formulae (PU]), ( P7D 
and ( pSD one obtains: A = sii? 26, B = {CP and T violation effects vanish, in the two 
flavor case) and C = cos 20 sin^ 20/3 in agreement with the general result. 
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Appendix C: Masses and mixing for large matter effects 



The dependence of the effective squared mass differences and mixing parameters in matter 
on the product pE^ is shown in figures ^, ^ and |^ (in the figures we assumed an electron 
fraction of |, therefore the density p and the potential V are simply proportional). Several 
important features are clearly visible: 

1. The squared mass eigenvalues and of the mixing parameters in matter are in general 
very different for neutrinos and anti- neutrinos. Reversing the sign of Amjs, the 
behaviour of v and I7's to a good approximation is simply exchanged. 

2. All four mixing parameters (the tree angles and the phase) change when matter is 
present, however, when the two mass scales Am^2 ^^^d |Am23| are of different orders 
of magnitude, as indicated by the the data, the angle 6*23 and the phase 5 remain 
approximately constant while the angles 612 and change much more dramatically. 

3. We can recognize several ranges of the parameter a = 2EyV where the behaviour of 
the solution has different characteristics: 

(i) When 2EyV -C Am^j; matter effects are negligible and the oscillations de- 
velop as in the vacuum case. 

(ii) For 2Eiy <^ |Am23| the effective mass and mixing of the state remain un- 
changed, but Am^2 6*12 can be modified by matter effects (this is interesting 
only if Am^2 ^ I ^"^23 1)- ^^^^ situation it is a a good approximation to use 
the well known two-flavor formulas to obtain (Amf2)™ and 9^ as a function of 
2VE^. 

(iv) When 2 E^V ~ |Am23| there is a resonance (6'^^ becomes |). The resonance is 
present for neutrinos if Am23 positive or for anti-neutrinos if Am|3 is negative. 

(v) When 2VE^ ^ |Am23|, the behaviour of the mixing angles and squared mass 
eigenvalues takes a simple form that will be discussed below. 

In the study of the effective parameters for large matter effects {\dim2^\/ {2V Ey) 0) 
one has to distinguish two cases: 

• Case A corresponds to neutrinos for Am23 positive or to anti-neutrinos for Am23 
negative. In this case one has: 

(Am?2)- Am^g, 



(Am^g)"^ 2VEy, 



2 



2 



sin2 ^ 1, cos^ ^-3 ^ sin^ ^13 ' ^'"^s 



2E,y 



sin^ 0^2 1, cos^ ff^^ const ~ sin' ^12 

A?Ti23 

Thin means that becomes asymptotically a massive pure state, therefore 
6^13 I and cos' ^13 vanishes rapidly (oc 2i?i^l^)~'. The angle 6^12 tends to a constant 
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value, reflecting the fact that the flavor is "sucked away" at the same rate from 
the \v\) and \v2) states. The squared mass difference (Am^a)™ grows without bound 



lEyV^ while (Am^g)™ approaches the constant value (— Am^ 



23/- 

Case B corresponds to neutrinos for Am^a negative or to anti-neutrinos for lS.m\^ 
positive. In this case one has: 

(Am^g)- ^ Am^3, 
sin2^-3-.sin2^^i3 (^)'' cos^^S-l, 



sin^r,-sin^^^i,('^) , cos^^- 



In this case it is |z/™) that becomes the massive pue state, therefore both and Q\2 
asymptotically vanish oc {2EyV^~^ ^ and it is (Am^g)™ ^^at grows large (— lE^V^ while 
Am23 remains approximately unchanged. In both cases the modifications to the the angle 
^23 and the phase b are small and vanish for Amf2/I^'^23l small. 

Note that the behaviour of the effective masses and mixing parameters is strikingly 
different for z/'s and z7's. However there are some important cancellations. For example 
collecting the expressions for the different factors one can verify that the the effective 
Jarlskog parameter in matter (Jm = (c^^)^ s^^ s^^ cj^ sin^m) is equal for v and V: 

•^ra ~ '^ra ^ -^vacuum ^^j^ y^^2 ' V 

(see for numerical calculation in fig. |^ as it is the case for the product of the three effective 
squared mass differences: 

{^ml,r^ {Aml,Y^ {Aml,r^ = (Am?^)^ i^rnl,)l (Aml.Z - Am^3 i2VE^y (92) 

Combining the two results one satisfies equation (|45|). 
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Figure 1: Electron neutrino fluence {d(j)u,J dE^) in a neutrino factory machine. The different curves are 
calculated for a fixed number of unpolarized /x decays, with different energy: = 5, 10, 20 and 40 GeV. 
Note that increasing E^ the integrated u fluence increases oc E'^, but the fluence at low energy decreases. 
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Figure 2: Electron neutrino fluence in a neutrino factory machine {d(j)v^ / dlog E^) . 
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Figure 3: Graphical representation of the flavor components of the neutrino mass eigenstates. The mixing 
angles have the same values in all three panels: 612 = 023 = 45°, ^13 = 12°, the phase S is 0, 180° and 
±90° in the left, center and right panel. Note how the quantities |(f/i,T|z^i,2)P depend on the value of \6\. 




Figure 4: Geometrical relation between the flavor and mass eigenvectors in the case of a real mixing 
matrix {6 = or n). The figure shows the projections in the (^'^,^'^-) plane of the mass eigenvectors. The 
fe vector is coming out of the plane of the figure. The parameters of the mixing matrix are indicated in 
the plots. See text for more discussion. 
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Figure 5: Effective squared mass eigenvalues in matter, plotted as a function of p, The top (bottom) 
panel is for i/ (F). The squared mass values are: = 0, rn| = 7 x 10~^ eV^, m| = 3 x 10~^ eV^; the 
mixing parameters are: 612 = 40°, ^23 = 45°, O13 = 7° , 6 = 45°. In this figure the elkec 
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Ej, p (GeV g cm ^) 

Figure 6: Mixing parameters in matter as a function of the product p: top panel sin^ "^Of^- bottom 
panel sin^26'5^. The squared mass values are: m\ = m\ = 7 x 10~^ eV^, m§ = 3 x 10~^ eV^; the 
mixing parameters in vacuum are: 6\2 = 40°, 623 = 45°, ^13 = 7°, ^ = 45°. The solid curves are 
for i^'s, the dashed curves for 77's Note the simple asymptotic forms of the parameters for large Ei, p. 
The dotted hues in the top panel is sin^ 20i2 (Amf2/(2£'j/V^))^. The dotted line in the bottom panel is 
sin^ 26*13 (Am?3/(2i;^y))2. 
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Figure 7: Effective mixing parameters in matter plotted as a function of the product Ei, p. The top panel 
shows sin^ 0^, the bottom panel shows the Jarlskog parameter: J = (cj^)^ sfs sf^ c^^ sin (5m- The 

squared mass values are: rn^ = 0, rr^ = 7 x 10~^ eV^, ml = 3 x 10~^ eV^; the mixing parameters in 
vacuum are: 6*12 = 40°, Q2?, = 45°, 9 is = 7°, S — 45°. The solid curves are for i^'s, the dashed ciu'ves for 
F's The dotted line in the bottom panel is Jyac ^"^^is Amf2/(2-Ei/l^)^ and indicates the asymptotic form 
of the parameter for large E^, p. Note how the asymptotic form is the same for i/'s and F's. 
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Figure 8: Oscillation probability for the transition — > plotted as a function of E^, for a fixed value 
of the pathlength L = 730 Km. The oscillation parameters are fixed, and are given inside the figure. The 
four curves are for neutrinos and anti-neutrinos in vacuum and in matter. 
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Figure 9: Oscillation probability for the transition plotted as a function of Ei, for a fixed value 

of the pathlength L = 3000 Km. The oscillation parameters are fixed, and are indicated in the figure. 
The solid (dashed) curves are for v (V). The top (bottom) panel gives the probability for propagation in 
vacuum (matter). 
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Figure 10: Oscillation probability for the transition v^, plotted as a function of for a fixed value 

of the pathlength L = 7000 Km. The oscillation parameters are fixed, and are indicated in the figure. 
The solid (dashed) curves are for u iv). The top (bottom) panel gives the probability for propagation in 
vacuum (matter). 
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Figure 11: In this figures we plot as a function of E^, (for a fixed value L = 730 Km) the product 

P{iye — > t^fj.) X E^. The different curves correspond to different values of the phase S. The values of the 
other parameters is indicated inside the plot. 
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Figure 12: In this figure we plot as a function of the neutrino energy (for a fixed value L = 730 Km) 
the product P{ve x E^. The different curve describe the probability with and without matter 

effects. 
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Figure 13: In this figures we plot as a function of the neutrino energy E^, the product P{i'e '^ii) x E^, 
for a fixed value of the neutrino pathlength L = 3000 Km. The top panel is for vacuum oscillations the 
bottom one includes matter effects. 
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Figure 14: The top panel shows as a function of the v pathlength L, the energy iJpeak where the proba- 
bihties for Ve Vp,^T transitions have the matter enhanced maximum. The bottom panel show the size 
of the matter enhancement. The different curves correspond to differcmt values of S13. For small S13 the 
position of the peak and the size of the enhancement are independent from S13. 
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Figure 15: In this figure are indicated some interesting regions in the space {L,Ei,) of the oscillation 
probability for fe ^ and Ve <-> Vt transitions. The region above the dot-dashed line is the one where 
the high energy expansion is valid. In the region delimited by the dashed line the oscillation probabilities 
in matter and in vacuum are to a good approximation equal. The line labeled A shows the v energy 
where the oscillation probability has the largest matter induced enhancement. The line is drawn only if 
the maximum enhancement is larger than 20%. The line labeled with a shows the relation L = 2V~^ 
for the Earth's crust (p = 2.8 g cm~^). The lines b and c show the relations Ei^ = |Am23| -^/(^tt) and 
= Amfa -^/(27r), that is the highest energy where the oscillation probability has a maximum. Line d 
indicates the approximate energy for which 613 passes through a resonance, and line e shows the energy 
above which matter effects modify significantly the mixing parameters. 



46 



